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Discussion 


The  p-version  of  Che  finice  element  method  is  a  new  approach  to  finite 
element  analysis  which  has  been  demonstrated  to  lead  to  significant  computa¬ 
tional  savings,  often  by  orders  of  magnitude  (This  approach  was  formerly 
called  the  constraint  method;  the  new  term  p-version  is  more  descriptive). 
Conventional  approaches  (called  the  h-version)  generally  employ  low  order 
polynomials  as  basis  functions.  Accuracy  is  achieved  by  suitably  refining 

the  approximating  mesh.  The  p-verslon  uses  polynomials  of  arbitrary  order 

0 

p  >  2  for  problems  in  plane  elasticity  where  C  continuity  is  required  and 
polynomials  of  order  p  >_  5  for  problems  in  plate  bending  where  continuity 
is  required. 

Hierarchic  elements  which  Implement  the  p-version  efficiently  are  used 

together  with  precomputed  arrays  of  elemental  stiffness  matrices. 

A  major  result  was  obtained  in  the  previous  reporting  period  on  the 

convergence  of  the  p-version  of  the  finite  element  method  for  the  case  of 

Co  continuity:  let  N  be  the  number  of  degrees  of  freedom.  In  polynomial 

regions,  if  the  h-version  converges  with  order  Q(l/Na)  then  the  p-version 

2a 

converges  with  order  0(1/N  )  provided  that  all  vertices  of  the  polygon  are 

vertices  of  the  triangulation. 

The  theoretical  order  of  convergence  of  the  p-version  in  the  case  of 
C*  continuity  has  still  not  been  determined  mathematically,  although  empi¬ 
rical  results  indicate  that  it  is  much  more  rapid  than  that  of  h-convergence. 
The  basic  problem  of  p-convergence  in  the  C^-  case  was  Investigated  during 
this  reporting  period.  By  using  modified  Bernstein  polynomials  defined  over 
triangles,  it  was  shown  that  the  p-version  does  Indeed  converge  in  the  energy 
norm,  although  the  optimal  rate  of  convergence  is  unknown.  The  construction 
of  these  modified  Bernstein  polynomials,  of  Interest  in  its  own  right,  is 
presented  in  [1]. 
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In  the  case  of  finite  elements  with  curved  boundaries,  the  standard 
h-version  procedure  is  to  use  isoparametric  mappings  together  with  numerical 
quadrature  to  compute  the  stiffness  matrix.  In  the  p-version  large  elements 
are  used  and  numerical  quadrature  becomes  inefficient  and  expensive.  In¬ 
stead  specific  mappings  are  sought,  which  map  curved  triangles  into  the 
standard  triangle,  in  such  a  way  that  the  resulting  elemental  stiffness 
matrix  can  be  computed  using  closed  form  integrals  and  pre-computed  arrays. 
Such  mappings  have  been  developed  in  [2]  for  triangles  with  one  curved  side 
which  is  either  elliptical  or  parabolic.  The  computer  implementation  of 
these  elements  will  greatly  increase  the  capabilities  of  the  p-version  for 
stress  analysis. 

Development  of  a  hierarchic  family  of  arbitrary  polynomial  order  for 
a  square-pyramidal  element  is  now  in  progress. 
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Abatracta  of  Papers  Presented  at  Meetings 

2 

[1]  Smooth  Approximation  to  a  Function  io  HQ(D)  by  Modified  Bernstein 
Polynomials  Ovar  Triangles  by  A.  G.  Kassos  and  I.  N.  Katz,  presented 
at  tha  SIAM  Fall  1979  Meeting,  November  12-14,  in  Denver,  Colorado* 

2 

Let  f  be  defined  on  a  bounded  polygonal  domain  D  in  R  and  let  S^T^}, 
1"1,.*.,N  be  a  triangulation  of  D.  Using  triangular  coordinates,  a  Bern¬ 
stein  polynomial  approximation  to  f  of  degree  p  >_  0  can  be  constructed 
over  each  T1*  Denote  this  approximation  by  B(f,p,Ti).  It  la  easily  shown 
that  the  global  approximation  induced  by  BCf,?,^)  i-1, . . . ,N  is  continuous 
on  D.  It  is,  however,  not  smooth  on  D.  In  this  paper,  we  construct  a  mod¬ 
ified  Bernstein  polynomial  approximation  on  each  T^,  denoted  by  B*(f .p.T^) , 
for  p  >_  5,  when  feC^(D)  and  we  show  that  the  global  approximation  induced 

by  B^f.p.T^)  i-l,...,N  is  smooth.  This  result  can  be  used  to  prove  that 
2 

any  fcH^CD)  can  be  approximated  smoothly  by  a  function  which  is  a  piecewise 
polynomial  of  sufficiently  high  degree  p  over  a  fixed  trlangulatlon  S.  The 
procedure  in  which  S  la  fixed  while  p  is  increased  is  called  the  p-verslon 
of  the  finite  element  method,  and  its  computer  implementation  has  been  de¬ 
scribed  elsewhere. 

[2]  Triangles  with  One  Curved  Side  for  the  P-version  of  the  Finite  Element 
Method,  by  I.  N.  Katz,  presented  at  the  SIAM  Spring  1980  Meeting, 

June  5-7,  in  Alexandria  Virginia. 

In  the  p-verslon  of  the  finite  element  method,  the  triangulation  of  the 
domain  la  fixed  and  the  degree  p  of  the  approximating  polynomial  is  varied. 

In  this  approach  the  triangles  are  generally  large  and  numerical  quadrature, 
in  tha  case  of  curved  elements,  may  be  expensive.  In  order  to  circumvent 
this  difficulty  we  consider  the  triangles  with  one  curved  side  in  which  a 
particular  shape  is  specified  for  the  curved  side  (c.g.  circular,  parabolic, 
elliptic).  A  mapping  of  the  curved  element  into  a  standard  element  is  sought 
with  the  property  that  the  resulting  elemental  stiffness  matrix  can  be  in¬ 
tegrated  in  closed  form  (without  recourse  to  numerical  quadrature).  Such 
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nappings  arc  given  for  the  following  two  specif icationa  of  the  curved  and 
straight  sides:  (1)  curved  side  parabolic,  straight  sides  arbitrary  (2) 
curved  side  elliptic,  straight  sides  parallel  to  the  axes  of  the  ellipse. 
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